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A  TWO-WAVE  EMPIRICAL  MODEL  OF  THE  PRESSURE 
FIELD  IN  A  LIQUID-FILLED  CYLINDRICAL  SHELL 

1.  INTRODUCTION 

A  comparison  of  experimental  data  with  a  model  of  the  interior  pressure  field  of  a 
liquid-filled  cylindrical  shell  often  shows  disagreement  between  prediction  and  measurement. 
One  reason  is  that  the  models  for  this  field  [1-3]  typically  do  not  take  into  account  such 
unmodeled  processes  as  internal  hydrophone  mount  noise,  support  cable  interaction,  and  other 
noise  effects  that  cause  deviations  from  the  measurements.  In  this  report,  a  model  of  the 
pressure  field  in  a  liquid-filled  cylindrical  shell  is  developed  based  on  experimental  pressure 
measurements  of  the  field.  The  underlying  assumption  is  that  two  waves  (an  extensional  and  a 
breathing  wave)  are  propagating  in  both  axial  directions  of  the  shell.  A  four-term  model  that 
describes  this  motion  is  derived,  and  the  six  unknown  model  constants  are  determined  with 
two  inverse  methods.  The  model  makes  no  assumptions  about  the  boundary  conditions  of  the 
shell  and  the  material  properties  of  the  shell  or  liquid.  Concerning  the  shell/liquid  interaction, 
the  only  assumption  is  that  the  dynamics  of  this  process  creates  a  breathing  wave  in  the 
structure  that  originates  at  the  bulkheads.  This  empirical  model  is  developed  as  a  "bridge" 
between  theory  and  experiment,  and  it  will  identify  the  frequency  ranges  where  disagreement 
typically  occurs. 


2.  SYSTEM  MODEL 


The  pressure  field  inside  the  liquid-filled  shell  is  derived  from  two  wave  equations  in 
the  spatial  domain,  both  of  which  use  pressure  as  the  independent  variable.  The  first  equation 
models  the  extensional  wave  contribution  and  is  written  as 


d^PJx,(o)  .  .2 


dx" 


+  kgPg{x,(o)  =  Q 


(1) 
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where  Pg(x,co)  is  the  temporal  Fourier  transform  of  the  pressure  that  is  generated  by  the 
extensional  wave,  x  is  the  spatial  location  (m),  w  is  the  excitation  frequency  (rad/s),  and  is 
the  complex  extensional  wavenumber,  which  is  equal  to 


where  Cg  is  the  complex  extensional  wave  speed  (m/s).  The  second  wave  equation  models  the 
breathing  wave  contribution  and  is  written  as 

+  =  0  ,  (3) 

dx 

where  Pfj{x,0))  is  the  temporal  Fourier  transform  of  the  pressure  that  is  generated  by  the 
breathing  wave  and  kfj  is  the  complex  breathing  wavenumber  (rad/m),  which  is  equal  to 

kt=—,  (4) 

where  is  the  complex  breathing  wave  speed  (m/s).  The  solutions  to  equations  (1)  and  (3) 
are  complex  exponential  functions,  and  they  can  be  added  together  using  the  principle  of 
superposition,  which  yields  the  total  pressure  in  the  shell  as 

P(x,co)  =  Pbix,(0)  +  Peix,0))  =  +  B{(o)q~^  +  C(6})e^«^  +  D{(0)t~^  ,  (5) 

where  P(x,co)  is  the  temporal  Fourier  transform  of  the  pressure  that  is  generated  by  both  the 
extensional  and  breathing  waves;  i  is  the  square  root  of  -1;  and  A{co),  C(<o),  and 

D((o)  are  wave  propagation  coefficients  determined  by  the  boundary  conditions.  It  is  now 
noted  that  the  pressure  field  at  x  divided  by  the  forward  accelerometer  is 

where  U  is  the  temporal  Fourier  transform  of  the  acceleration  at  the  forward  end  of  the  shell; 
A(co),  B{co),  C(a)),  andD(6))  are  wave  propagation  coefficients;  and  P(x,co)IU  has  units  of 
Pa//(m/s^).  Physically,  the  coefficient  A(q})  corresponds  to  forward-traveling  breathing  wave 
energy,  B((o)  corresponds  to  aft-traveling  breathing  wave  energy,  C((0)  corresponds  to 
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forward- traveling  extensional  wave  energy,  and  D^co)  corresponds  to  aft- traveling  extensional 
wave  energy. 

There  are  six  unknowns  in  equation  (6),  which  are  the  four  wave  propagation 
coefficients  and  the  two  forced  (breathing  and  extensional)  wavenumbers.  Using  a  laboratory 
configuration  in  which  a  longitudinal  shaker  is  placed  at  the  forward  end  of  a  shell  allows  all 
these  unknowns  to  be  estimated.  Due  to  the  nonlinear  (and  sinusoidal)  nature  of  the  parameters 
in  equation  (6),  the  most  numerically  stable  method  is  to  first  estimate  the  extensional 
wavenumber,  then  the  breathing  wavenumber,  and  finally  the  wave  propagation  coefficients. 
A  method  to  estimate  the  extensional  and  breathing  wavenumbers  has  previously  been 
developed  [4,5],  and  a  shortened  version  is  included  here  as  section  3. 


3.  ESTIMATION  OF  THE  FORCED  WAVENUMBERS 
In  order  to  estimate  the  wave  propagation  coefficients  and  the  breathing  wavenumber 
(and  the  corresponding  breathing  wave  speed),  it  is  necessary  to  know  the  extensional 
wavenumber  (which  is  related  to  the  corresponding  extensional  wave  speed  by  equation  (2)). 
This  complex,  frequency-dependent  quantity  is  determined  using  measurements  from  forward 
and  aft  impedance  heads  that  are  attached  to  the  shell  ends  in  the  test  facility  [4].  Although 
these  measurements  contain  a  breathing  wave  contribution,  it  typically  occurs  only  at  low 
frequencies  and  can  be  easily  discempd  from  the  extensional  wave  effects.  Such  behavior  is 
described  in  the  experiment  section. 

The  governing  differential  equation  of  the  extensional  wave  is  expressed  on  the  spatial 
domain  as  a  single  wave  equation  with  particle  displacement  as  the  independent  variable: 


d^U{x,(o) 

dx2 


+  keUix,(O)  =  0  , 


(7) 


where  U(x,(o)  is  the  temporal  Fourier  transform  of  the  axial  displacement.  Equation  (7)  is  a 
one-dimensional  "lumped"  approximation  of  the  extensional  wave  motion  in  the  system  that 
combines  the  effect  of  the  shell  and  the  liquid  together  as  a  homogeneous  medium  that  supports 
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longitudinal  wave  motion.  Although  this  approximation  is  not  sufficient  to  model  the  radial 
motion  of  the  shell,  it  is  an  accurate  model  of  axial  motion  and  the  corresponding  extensional 
wave  propagation  in  the  liquid-filled  shell.  The  energy  attenuation  in  the  shell  is  defined  with  a 
structural  damping  law,  and  therefore  the  wave  speed  is  a  complex  quantity.  The  real  part  of 
the  wave  speed  corresponds  to  energy  transmission  and  the  imaginary  part  corresponds  to 
energy  attenuation. 

The  solution  to  equation  (7)  is 

U{x,(0)  =  ,  (8) 

where  G(g))  and  H{co)  are  coefficients  determined  by  the  boundary  conditions  at  the  ends  of  the 
shell.  The  temporal  Fourier  transform  of  the  axial  force  in  the  shell  is 

where  is  the  cross-sectional  area  of  the  shell  (m^)  and  is  the  effective  longitudinal 
modulus  of  the  shell  (N/m^).  The  known  parameters  in  equations  (8)  and  (9)  are  the  location 
of  the  sensors  (;c)  and  the  frequency  of  excitation  (co).  Although  the  effective  longitudinal 
modulus  (E^)  is  unknown,  equation  (9)  will  be  rewritten  as  a  ratio  of  forces  permitting  the 
cancellation  of  this  term  and  the  cross-sectional  area  (A^).  Additionally,  the  coefficients  G(co) 
and  Hico)  are  unknown;  however,  they  will  be  condensed  out  of  the  mathematical 
relationships.  The  inversion  of  equations  (8)  and  (9)  at  the  sensor  locations  will  allow  for  a 
measurement  of  the  unknown  extensional  wavenumber  and  extensional  wave  speed  Cg.  This 
technique  is  described  next. 

In  the  extensional  wave  measurement  part  of  the  experiment,  the  forward  and  aft  pair  of 
sensors  (impedance  heads)  described  earlier  collect  data  that  are  in  the  form  of  transfer 
functions  between  each  pair.  The  position  of  the  forward  pair  of  sensors  is  defined  as  jc  =  0 
and  of  the  aft  pair  as  x  =  L,  where  L  is  the  length  of  the  shell  (m).  The  two  transfer  function 
measurements  used  are  the  forward  displacement  divided  by  the  aft  displacement  and  the 
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forward  force  divided  by  the  aft  force.  Their  theoretical  form  can  be  rewritten  using  equations 


(8)  and  (9)  as 


G{(o)  +  H{co) 


U(L,  CO)  C(G))e“'e^  +  H(co)c~'^^ 


-ik.L  ~ 


F(0,co)^  G{C0)-H{(0) 

FiL,(o)  G(ty)e^«^-//(6))e"'^«^ 


where  and  R2  are  transfer  function  data  from  the  experiment  at  a  specific  test  frequency. 
Equations  (10)  and  (11)  are  rewritten  as  functions  of  H(co)/G(cd)  and  are  set  equal  to  each 
other,  yielding 

cos(k,L)  =  ^^  =  (f>  ,  (12) 

R2+R1 

where  ^is  a  complex  quantity.  Using  an  angle-sum  relationship  on  the  complex  cosine  term  in 
equation  (12)  and  separating  the  equation  into  real  and  imaginary  parts  results  in 


cosh[Im(^g)L]  = 


Re(0) 

cos[Re(^g)L] 


sinh[Im(/:g)L]  =  — 


-lm(^) 


"  sin[Re(kg)L] 

where  Re  denotes  the  real  part  and  Im  denotes  the  imaginary  part  of  the  corresponding  complex 


quantity. 


Equation  (14)  is  now  squared  and  subtracted  from  the  square  of  equation  (13),  yielding 


{cosh[Im(^g)L]}^ - {sinh[Im(/:g)L]}^  =  - — ^ ^  ^ ^  =  1  .  (15) 

{cos[Re(A:g)L]}^  {sin[Re(/t^)L]} 

Equation  (15)  can  be  simplified  with  trigonometric  power  relationships  to 
cos[2Re(/:g)L]  = 

[Re(0)f  -h -  J([Re(<p)f  +  [lm(0)]^)  - (2[Re(0)]^  - 2[Im((p)f  -l)=s.  (16) 
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Note  that  only  a  negative  sign  in  front  of  the  radical  is  used.  The  real  part  of  in  equation 
(16)  is  now  solved  for  by 

— Arc  cos(s)  +  —  n  even 
2L  ^  2L 

Re(ite)  =  -^  ,  (17) 

1  .  ,  .  nn  ,, 

— Arc  cos{-s)  +  —  n  odd 

[21  2L 

where  n  is  a  nonnegative  integer  and  capital  A  denotes  the  principal  value  of  the  inverse  cosine 
function.  The  value  of  n  is  determined  from  the  function  s,  which  is  a  cosine  function  with 
respect  to  frequency.  At  zero  frequency,  n  is  0.  Every  time  s  cycles  through  it  radians,  n  is 
increased  by  1.  The  imaginary  part  of  kg  is  determined  by  adding  equations  (13)  and  (14) 
together,  resulting  in 


Im(/:e)  =  —  loge 


Re(^) 

cos[Re(^g)L] 


lm(0) 

sin[Re(^g)L] 


Now  that  the  real  and  imaginary  parts  of  the  wavenumber  kg  are  known,  the  complex-valued 
extensional  wave  speed  can  be  determined  at  each  frequency  by 


Cg  =  Re(Cg)  -f-  ilm(Ce)  =  ^  .  (19) 

ke 

Note  that  the  extensional  wave  speed  has  been  measured  without  knowing  the  boundary 
conditions  at  x  =  0  and  x  =  L. 

Once  the  extensional  wave  speed  is  known,  five  independent  equally  spaced 
measurements  of  the  spatial  pressure  field  are  needed  to  eliminate  the  wave  propagation 
coefficients  and  solve  for  the  breathing  wave  speed  [5].  Without  loss  of  generality,  the  origin 
of  the  coordinate  system  is  defined  as  x  =  0  at  the  middle  (third)  pressure  sensor  (hydrophone). 
Equation  (6)  is  written  to  correspond  to  the  locations  of  the  five  pressure  sensors  as 

^  b{co)q^^'^  +  C{(o)Q~^^  -t- D(£0)e^*^^  =  Si  ,  (20) 

=  A{o))c~^  +  5(<»)e^«^  +  C{0))Q~^  +  D(co)e'**^  =  Si  ,  (21) 


6 


TR  11,143 


P(0,O)) 


=  A(co)  +  B((o)  +  C{(o)  +  D{(0)  =  53  , 


+  B{(o)q~^^  +  C((W)e^*^^  +  D{,(o)t  ,  (24) 

where  8  is  the  sensor-to-sensor  spacing  (m)  and  Sj  through  correspond  to  the  measured 
transfer  function  data  of  the  hydrophone  divided  by  the  forward  accelerometer  at  a  specific 
frequency. 

Equations  (21)  and  (23)  are  now  added  together  to  yield 

(A  +  B)cos(ke5)  +  (C+D)cosikt,d)  =  (l/2)(52  +  54)  ,  (25) 

and  equations  (20)  and  (24)  are  added  together  to  produce 

(A  +  B)cos(ke25)  +  (C  +  D)cos(ki,2S)  =  (l/2)(5i  +  S5)  .  (26) 

Equation  (22)  is  then  rewritten  with  the  term  (A  +  B)  on  the  left-hand  side  and  is  substituted 

into  equations  (25)  and  (26),  yielding 

(l/2)(52-H54)-53Cos(^,g) 
cosikijS)  -  cos(kg5) 


cos(«^2d)  -  cos(kg2d) 

respectively.  Equations  (27)  and  (28)  are  next  set  equal  to  each  other.  Applying  a  double 
angle  trigonometric  relationship  to  the  cos(^/,25)  term  then  produces 


X  cos^  (kjjS)  +  Y  cos(ki,S)  +  Z  =  0  , 


where 


X  =  (52  S4 )  ^  2S3  cosikgS)  , 

F  =  53  cosikgS)  -  (1  /  2)(5i  -I-  55 )  , 


(29b) 
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and 


Z  =  [(l/2)(5i  +  Ss)  +  S2]cos(ke5) - (l/2)(52  +  54)008(^^25) -  (l/2)(52  +  -^4)  •  (29d) 

Equation  (29)  is  a  quadratic  form  with  the  solution 

,  (30) 

where  y^is  a  complex  quantity.  Rewriting  equation  (30)  as  real  and  imaginary  terms  produces 
cos[2Re()t^)5]  = 


[Re(v^)]2+[Im(vA)]2 


-  ^([Re(i/)f  +  -  (2[Re(yr)f  -  2[Im(i/^)f  - 1)  =  r  ,  (31) 


where  Re  denotes  the  real  part  and  Im  denotes  the  imaginary  part  of  the  corresponding  complex 
quantity.  Note  that  only  a  negative  sign  in  front  of  the  radical  is  used  in  equation  (31). 
However,  both  the  negative  and  the  positive  signs  in  front  of  the  radical  in  equation  (30)  are 
needed.  The  real  part  of  kj^  in  equation  (3 1)  is  now  solved  for  by 

1  .  .  .  mn 

— Arc  cos(r)  +  — —  m  even 

2o  2o 


Re(4)  = 


1^— Arccos(-r) 


+ 


mn 

Id 


m  odd 


(32) 


where  m  is  a  nonnegative  integer  and  capital  A  denotes  the  principal  value  of  the  inverse  cosine 
function.  The  value  of  m  is  determined  from  the  function  r,  which  is  a  cosine  function  with 
respect  to  frequency.  At  zero  frequency,  m  is  0.  Every  time  r  cycles  through  n  radians,  m  is 
increased  by  1.  The  imaginary  part  of  is  determined  from  equation  (30),  resulting  in 


Im(^^,)  = -jlogej 


Re(v/') 


cos[Re(/:^)5] 


Im(v^)  1 

sin[Re(^^)5]J 


(33) 


Now  that  the  real  and  imaginary  parts  of  the  wavenumber  kfj  are  known,  the  complex-valued 
breathing  wave  speed  can  be  determined  at  each  frequency  with 


(O 


Ci,  =  Re(c^)  +  iIm(c^)  =  — 

4 


(34) 
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Use  of  this  method  produces  two  wave  speed  measurements  because  of  the  retention  of  the 
positive  and  negative  signs  in  equation  (30).  One  of  the  wave  speeds  is  the  breathing  wave 
speed,  and  the  other  is  the  extensional  wave  speed,  which  was  previously  known.  The 
extensional  wave  speed  is  typically  at  least  one  order  of  magnitude  greater  than  the  breathing 
wave  speed. 

4.  ESTIMATION  OF  THE  WAVE  PROPAGATION  COEFFICIENTS 

Two  methods  are  now  used  to  estimate  the  wave  propagation  coefficients.  The  first, 
which  uses  four  hydrophone  measurements  and  (mathematically)  produces  an  exact  solution 
for  the  coefficients,  is  called  the  direct  inverse  method.  The  second  uses  N  hydrophone 
measurements  and  a  Gauss  linearization  scheme  with  an  ordinary  least  squares  constraint  to 
solve  for  the  coefficients  [6].  Both  these  methods  are  mathematically  developed  in  this  section 
and  applied  to  data  in  the  experiment  section. 

In  the  first  method,  equation  (6)  is  rewritten  at  four  different  hydrophone  measurement 
locations  as 

=  A(a))e^^^i  +  fi(<u)e"'^*^l  +  C(cy)e^«"'i  +  D(Q))e~^  =  Ti  ,  (35) 

=  A(co)c^b^^  +  B((o)t-^bX2  +  c(a))e^^^2  +  £)((u)e-^^^2  =  ,  (36) 

=  A(a))e^^^3  +  +  c(a))e^«^3  +  =  Tj  ,  (37) 

and 

=  A(G))e'*^^4  +  +  C(a))e^«^4  +  =  T4  ,  (38) 

where  Tj,  T2,  T^,  and  T4  are  the  transfer  function  data  at  a  specific  frequency.  Unlike  the 
estimation  of  the  breathing  wave  speed,  it  is  not  necessary  that  the  spatial  distances  JCi  through 
X4  be  equally  spaced.  Equations  (35)-(38)  are  now  written  in  matrix  form  as 
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Mxi 

Q-ikbxl 

e^exi 

Q-ikeXi  - 

'A(CO) 

ikbX2 

Q^bX2 

t^eX2 

Q-if^eX2 

B(0)) 

Q-ikbXi 

C^eX3 

Q-ikeX^ 

Cico) 

Q'i^bXA 

e^eXA 

Q^eXA 

Pico) 

The  coefficients  are  then  solved  for  by  multiplication  of  each  side  of  equation  (39)  by  a  matrix 
inverse,  which  results  in 

iTi  ■ 

B{(0)  e^bX2  Q-'^bX2  Q^e^2  Q^eX2  Ty 

^  =  i  i.  ^  (40) 

Cico)  t^bx^ 

p{(0)]  Q^bH  Q-^bH  e'*e-'4  Q-^eXA  [T^ 

In  the  second  method,  or  the  Gauss  linearization  scheme,  the  matrix  iteration  equation 
to  estimate  the  wave  propagation  coefficients  is 


b(‘+»  =b<«+[[X<‘)f  [X»>]]''[[X<«f(y-h<«)]  , 


where  the  superscript  k  represents  the  /:th  iteration;  the  vector  b  is  the  wave  propagation 
coefficient  vector  and  is  equal  to 

'A(co)) 

B(co) 

b  =  i  (42) 

C(a)) 

P(Q)) 

X  is  the  following  by  4  sensitivity  coefficient  matrix 


Q^bxi 

q-'Hxi 

e^exi 

Q-^eXl 

Q^bX2 

Q-^bX2 

Q^eX2 

^■ikeX2 

^^bXN 

Q-^bXN 

Q^eXN 

Q-^eXN 

y  is  the  N  dimensional  vector  containing  the  measurement  data  as 
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y  = 


'Ti' 

Tl 


(44) 


and  h  is  the  iV  dimensional  vector  containing  the  model  calculated  using  the  current  value  of  the 


estimated  parameters  and  is  equal  to 


A((i»e^^^2  +c(tu)e^^^2  +D(to)e-^^^2 


(45) 


[A{(0)Q^bXN  +BiQ))c-^bXN  +C{co)c^eXN  +D(G»e-^^^^  J 

Once  the  equations  are  assembled,  the  matrix  b  can  be  determined  by  iteration  on  the  left-hand 


side  of  equation  (41)  at  each  specific  frequency. 


5.  EXPERIMENT 

Use  of  this  model  with  the  inverse  method  corresponds  to  the  physical  testing 
configuration  in  the  Axial  Vibration  Test  Facility  (AVTF)  at  the  Naval  Undersea  Warfare 
Center/Division  Newport,  as  shown  in  figure  1.  The  AVTF  has  been  designed  to  provide  a 
simple  procedure  for  testing  long  structures  under  varying  tensions  and  temperatures.  The 
longitudinal  shaker  at  its  forward  end  provides  axial  excitation  to  the  structure.  A  rope  attached 
to  the  aft  end  of  the  structure  and  a  winch  allows  the  tension  to  be  adjusted.  A  mass  is  attached 
between  the  shell  and  the  rope  to  increase  the  force  levels  and  decrease  the  acceleration  levels. 
This  mass  also  produces  an  impedance  change  at  the  end  of  the  bar  that  is  sufficiently  large  to 
allow  accurate  modeling  of  the  rope  behavior  by  a  spring  and  damper  rather  than  by  a 
continuous  media  expression.  A  rail  from  which  thin  Kevlar  lines  can  be  hung  to  provide 
lateral  support  to  heavy  or  long  test  specimens  runs  the  entire  length  of  the  facility.  The  unit  is 
completely  surrounded  by  an  air-conditioned  PVC  duct  to  permit  temperature-dependent 
testing.  To  collect  data  during  a  test,  impedance  heads  are  attached  to  the  forward  and  aft  ends 
of  the  structure  (as  described  earlier).  Each  impedance  head  consists  of  a  single  axial  force 
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transducer  and  an  accelerometer.  Additional  data  are  collected  by  hydrophones  in  the  liquid- 
filled  shell.  A  load  cell  that  measures  the  tension  on  the  structure  is  located  between  the  rope 
and  the  mass. 

Equipment  used  in  this  experiment  includes  a  Zonic  Corporation  model  1215-10-T- 
ZSP86  hydraulic  shaker,  two  PCB  Piezotronics  model  348A  accelerometers,  two  PCB 
Piezotronics  model  233A  force  transducers,  a  Lebow  model  31/3  load  cell,  and  a  Cordem 
Corporation  model  1215-RMO  tension  drum. 

In  the  first  part  of  the  experiment,  which  measures  the  extensional  wave  speed,  a 
longitudinally  reinforced  liquid-filled  urethane  shell  containing  five  equally  spaced 
hydrophones  was  placed  in  tension,  as  shown  in  figure  1.  The  shell  had  a  mean  radius  of 
0.015  m  and  a  thickness  of  0.0028  m;  the  internal  liquid  had  a  density  of  760  kg/m^.  The  axial 
tension  on  the  shell  was  890  N  and  the  stressed  length  was  12.0  m.  The  point  mass  had  a 
weight  of  13.6  kg.  The  data  from  the  force  transducers,  accelerometers,  and  hydrophones 
were  acquired  in  the  time  domain  with  a  Hewlett  Packard  (HP)  3562  dynamic  signal  analyzer. 
The  analyzer  then  Fourier  transformed  the  raw  data  to  the  frequency  domain  to  obtain  the 
desired  transfer  functions.  The  test  was  run  with  a  frequency  range  between  3  and  100  Hz. 

Equations  (7)-(19)  were  applied  to  the  experimental  test  data  taken  with  the  force 
transducers  and  accelerometers,  and  the  resulting  extensional  wave  speed  of  the  structure  was 
found.  Figure  2  shows  the  calculated  extensional  wave  speed  versus  frequency.  The  solid  line 
corresponds  to  the  computed  wave  speed  and  the  dashed  line  to  the  ordinary  least  squares 
(OLS)  straight-line  fit.  The  OLS  fit  was  applied  to  the  data  between  35  and  100  Hz  to 
minimize  the  effect  of  the  breathing  wave  interaction  seen  at  lower  frequencies.  The  resulting 
OLS  fit  was  Cg  =  677.3  +  0.8/  (m/s)  for  the  real  part  and  =  83.3  -  0.05/  (m/s)  for  the 
imaginary  part  (where/is  frequency  in  Hertz).  These  same  values  are  used  for  the  extensional 
wave  speed  in  the  following  calculations  in  order  to  determine  the  breathing  wave  speed  and 
estimate  the  wave  propagation  coefficients. 
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In  the  second  part  of  the  experiment,  transfer  functions  of  the  hydrophone  data  divided 
by  the  forward  accelerometer  data  were  collected.  The  five  hydrophones  were  spaced  at 
intervals  of  1.83  m  with  a  distance  of  2.96  m  from  the  forward  end  of  the  shell  to  the  first 
hydrophone.  Equations  (20)-(34)  were  applied  to  the  experimental  test  data,  and  the  resulting 
breathing  wave  speed  of  the  structure  was  calculated.  Figure  3  is  a  plot  of  the  breathing  wave 
speed  versus  frequency.  The  plus  symbols  mark  the  data  from  this  experiment  that  were 
obtained  using  the  method  developed  in  section  3.  In  the  figure,  it  is  noted  that  the  breathing 
wave  for  this  specific  structure  is  spatially  coherent  from  3  to  13  Hz.  Above  13  Hz,  the 
imaginary  part  of  the  measurement  begins  to  diverge.  It  is  also  possible  that  effects  other  than 
spatial  incoherence  are  preventing  the  breathing  wave  from  being  identified  at  higher 
frequencies.  Changing  the  extensional  wave  speed  by  ±20  percent  and  recalculating  the 
breathing  wave  speed  produced  a  change  in  the  breathing  wave  speed  of  less  than  1  percent. 
Thus,  the  breathing  wave  speed  measurement  method  is  very  insensitive  to  incorrect 
extensional  wave  speeds.  Additionally,  the  breathing  wave  speed  was  also  insensitive  to 
varying  axial  tensions. 
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Figure  2.  Extensional  Wave  Speed  Versus  Frequency 
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Once  the  extensional  and  breathing  wave  speeds  are  measured,  equations  (35)-(45)  can 
be  applied  to  the  hydrophone  data.  Although  four  hydrophone  measurements  are  sufficient  to 
define  the  empirical  pressure  field,  they  are  insufficient  to  verify  the  accuracy  of  the  model. 
This  is  because  the  model  (by  definition)  is  equal  to  the  four  data  points  used  to  create  it.  To 
resolve  the  problem  of  determining  the  accuracy  of  the  model,  12  hydrophones  were  placed  in 
the  shell  and  were  spaced  at  intervals  of  0.5  m  with  a  distance  of  3.0  m  from  the  forward  end 
of  the  shell  to  the  first  hydrophone.  The  coherence  versus  frequency  plots  for  each 
hydrophone  measurement  are  included  in  appendix  A.  The  signal-to-noise  ratio  of  each 
hydrophone  and  the  forward  accelerometer  were  approximately  60  dB.  The  models  were  fit 
using  two  different  variations  of  equation  (6):  the  first  included  the  A(o})  term  and  is  referred 
to  as  the  four-term  model  and  the  second  omitted  the  A(©)  term  and  is  referred  to  as  the  three- 
term  model.  This  A(o})  term  corresponds  to  the  breathing  wave  pressure  generated  by  the  aft 
bulkhead  of  the  shell.  The  effect  of  this  model  variation  is  discussed  below. 

Figures  4, 5,  6,  and  7  are  plots  of  the  data  and  two  models  of  the  pressure  field  versus 
distance  at  5.4,  12.7,  32.1,  and  61.2  Hz,  respectively.  The  pressure  field  versus  distance  at 
41  different  frequencies  is  included  in  appendix  B.  The  (pressure  field)  models  were  created 
using  the  direct  inverse  method  (equations  (35)-(40)).  The  solid  line  represents  the  pressure 
field  including  the  A(q))  term  and  is  based  on  the  spatial  locations  corresponding  to  the  first  (xj 
=  3.0  m),  fourth  (X2  =  4.5  m),  eighth  (x^  =  6.5  m),  and  twelfth  {x^  =  8.5  m)  hydrophone 
sensors.  The  dashed  line  shows  the  pressure  field  excluding  the  A((o)  term  and  is  based  on  the 
spatial  locations  corresponding  to  the  first  ix2  =  3.0  m),  sixth  (JC3  =  5.5  m),  and  twelfth  (x^  = 
8.5  m)  hydrophone  sensors.  The  x■^^  location  is  not  used  for  the  three-term  model  because  it 
corresponds  to  the  A(a))  term.  The  extensional  wave  speeds  used  in  deriving  both  models 
were  determined  from  figure  2,  and  the  breathing  wave  speeds  from  figure  3.  For  the  region 
above  15  Hz,  the  breathing  wave  speed  was  extrapolated. 
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Figure  4.  Pressure  Field  Versus  Distance  at  5.4  Hz  With 
Model  Determined  Using  the  Direct  Inverse  Method 
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Figure  6.  Pressure  Field  Versus  Distance  at  32.1  Hz  With 
Model  Determined  Using  the  Direct  Inverse  Method 
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Figures  8,  9,  10,  and  1 1  are  plots  of  the  data  and  two  models  of  the  pressure  field 
versus  distance  at  5.4, 12.7,  32.1,  and  61.2  Hz,  respectively.  The  models  were  created  using 
the  Gauss  linearization  method  (equations  (41)-(45)).  The  solid  line  shows  the  pressure  field 
including  the  A(ci))  term  and  is  based  on  all  12  hydrophone  measurements.  The  dashed  line 
depicts  the  pressure  field  excluding  the  A(co)  term  and  is  also  based  on  all  12  hydrophone 
measurements.  The  extensional  wave  speeds  used  in  creating  the  model  were  determined  from 
figure  2,  and  the  breathing  wave  speeds  from  figure  3.  For  the  region  above  15  Hz,  the 
breathing  wave  speed  was  extrapolated.  Due  to  the  similarity  of  the  Gauss  linearization  method 
to  the  direct  inverse  method,  no  additional  plots  of  the  Gauss  linearization  method  are  included. 

The  four  different  variations  of  the  model  were  all  computed  and  compared  to  the  data. 
The  formula  used  for  this  comparison,  in  decibels,  was 


where 


fiico)  =  20  logio 


f-L'ly 

hi(coi-\yiico) 

y(Q)) 

(46) 


1  ^ 

The  term  /  is  the  i'th  member  of  the  corresponding  vector  and  N  is  the  total  number  of 
measurements  {N  =  12).  The  average  difference  between  the  models  and  the  measurements 
(equation  (46))  at  various  frequencies  is  listed  in  table  1. 
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Figure  8.  Pressure  Field  Versus  Distance  at  5.4  Hz  With 
Model  Determined  Using  the  Gauss  Linearization  Method 
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Figure  9.  Pressure  Field  Versus  Distance  at  12.7  Hz  With 
Model  Determined  Using  the  Gauss  Linearization  Method 
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Table  1.  Average  Difference  Between  Model  and  Measurement  in  Decibels 


Frequency 

Direct  Inverse  Method 

Gauss  Linearization  Method 

(Hz) 

4  Term 

3  Term 

4  Term 

3  Term 

3.0 

0.6 

0.8 

0.6 

1.0 

5.4 

0.5 

0.8 

0.5 

2.3 

7.9 

0.4 

1.1 

0.4 

1.4 

10.3 

0.2 

0.7 

0.3 

2.4 

12.7 

0.6 

1.4 

0.4 

1.8 

15.1 

0.5 

1.2 

3.3 

0.7 

17.5 

1.2 

0.9 

3.2 

1.5 

20.0 

1.0 

1.4 

2.1 

2.5 

22.4 

0.9 

2.5 

3.7 

2.1 

24.8 

0.9 

4.2 

2.3 

2.5 

27.2 

1.2 

3.7 

2.0 

2.4 

29.7 

2.2 

3.2 

3.1 

3.5 

32.1 

3.6 

3.0 

4.8 

4.4 

34.5 

4.0 

2.6 

6.3 

3.9 

37.0 

5.1 

2.5 

4.5 

3.0 

39.4 

6.9 

2.5 

3.5 

2.7 

41.8 

4.8 

2.1 

2.8 

2.1 

44.2 

3.9 

1.8 

2.2 

2.0 

46.7 

3.3 

1.7 

2.1 

2.5 

49.1 

3.1 

1.8 

2.4 

2.6 

51.5 

2.9 

2.0 

2.4 

2.6 

53.9 

2.7 

2.1 

2.4 

2.7 

56.3 

2.7 

2.3 

2.5 

2.8 

58.8 

2.2 

2.0 

1.9 

2.2 

61.2 

2.3 

1.7 

1.7 

1.7 

63.6 

2.5 

1.8 

1.9 

1.8 

66.1 

2.4 

1.9 

2.0 

1.7 

68.5 

2.3 

1.9 

1.8 

1.7 
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Table  1.  Average  Difference  Between  Model  and  Measurement  in  Decibels  (Cont'd) 


Frequency 

Direct  Inverse  Method 

Gauss  Linearization  Method 

(Hz) 

4  Term 

3  Term 

4  Term 

3  Term 

70.9 

2.0 

1.9 

1.7 

1.7 

73.3 

1.9 

1.9 

1.8 

2.0 

75.8 

1.9 

1.6 

1.6 

1.6 

78.2 

1.6 

1.5 

1.3 

1.4 

80.6 

1.7 

1.1 

1.3 

1.3 

83.0 

1.3 

1.5 

1.4 

1.5 

85.4 

1.4 

2.2 

1.9 

2.1 

87.9 

1.5 

3.0 

2.2 

2.4 

90.3 

1.8 

4.1 

2.5 

2.6 

92.7 

1.7 

3.4 

2.5 

2.5 

95.2 

2.2 

3.2 

2.6 

2.6 

97.6 

2.1 

3.3 

2.6 

2.7 

100.0 

1.8 

3.2 

2.6 

2.7 

The  four-term  models  are  extremely  accurate  from  3  to  approximately  13  Hz.  At  13 
Hz,  the  Gauss  linearization  method  begins  to  diverge  from  the  direct  inverse  method.  The 
direct  inverse  method  is  still  accurate  up  to  approximately  27  Hz.  From  27  to  50  Hz,  there  is  a 
difference  of  about  4  dB  between  the  experiment  and  both  models.  Because  the  models  are 
being  defined  by  the  experimental  values  and  the  coherence  of  the  measurements  are  almost 
unity,  this  difference  is  larger  than  expected.  There  could  be  a  number  of  reasons  why  this 
divergence  occurs: 

•  A  transverse  wave  is  introducing  unmodeled  energy  into  the  measurements. 

•  The  wave  speeds  are  not  consistent  across  the  length  of  the  module. 

•  The  hydrophones  themselves  are  creating  significant  breathing  wave  energy  sources 
in  this  frequency  range. 
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•  The  transverse  support  lines  are  a  noise  source. 

•  The  spatial  coherence  length  of  the  breathing  wave  generated  at  the  aft  bulkhead  is 
causing  a  situation  where  some  breathing  wave  energy  is  being  input  into  the  aft  hydrophones 
but  not  into  the  forward  ones,  which  creates  an  inconsistent  model. 

•  Nonaxisymmetric  effects  of  wave  motion  are  being  detected  by  the  hydrophones. 

•  The  hydrophone  mounts  are  transmitting  unmodeled  energy  into  the  hydrophones. 

The  models  above  50  Hz  have  a  difference  from  the  measured  values  of  about  2  dB.  This  is 
considered  borderline  for  accuracy  in  this  particular  experiment. 

In  the  region  above  50  Hz,  the  four-term  models  may  be  more  accurate;  however,  they 
are  creating  pressure  levels  that  are  too  high  in  the  spatial  region  around  the  aft  bulkhead. 
Additionally,  the  phase  angle  associated  with  these  models  at  the  aft  bulkhead  is  incorrect  in 
this  frequency  range,  which  is  why  the  three-term  models  are  included  in  the  analysis. 
Because  the  three-term  models  assume  only  a  single-direction  breathing  wave  moving  from  the 
forward  bulkhead  to  the  aft  end  of  the  shell,  they  correspond  more  closely  to  the  physics  of  the 
problem  in  the  higher  frequency  range,  even  though  the  average  difference  between 
measurement  and  model  may  be  greater  than  that  of  the  four-term  models. 

6.  CONCLUSIONS  AND  RECOMMENDATIONS 

A  two-wave  empirical  model  of  the  pressure  field  in  a  liquid-fiUed  cylindrical  shell  can 
be  derived  based  on  four  hydrophone  measurements  in  the  liquid.  Two  comparably  accurate 
methods  were  used  to  derive  this  pressure  field:  a  direct  inverse  method  and  a  Gauss 
linearization  method.  A  three-  and  a  four-term  version  of  both  methods  was  investigated.  It 
was  found  that  at  the  frequency  range  of  3  to  30  Hz,  the  four-term  model  was  more  accurate 
and  in  the  frequency  range  of  50  to  100  Hz,  the  three-term  model  was  more  accurate.  It  is 
likely  that  an  unmodeled  noise  source  is  contributing  to  the  pressure  field  from  30  to  50  Hz. 


29 


TR  11,143 


The  experiment  should  be  rerun  using  shells  of  different  diameter,  stiffness,  and  length  to 
determine  if  these  model-to-experiment  differences  are  shell  dependent. 
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APPENDIX  A  -  COHERENCE  VERSUS  FREQUENCY  PLOTS 
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Figure  A-1.  Coherence  Versus  Frequency  for  Channel  1 
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APPENDIX  B  -  PRESSURE  FIELD  VERSUS  DISTANCE  PLOTS 


B-1 


Magnitude 


TR  11,143 


Figure  B-2.  Pressure  Field  Versus  Distance  at  5.4  Hz 
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Figure  B-4.  Pressure  Field  Versus  Distance  at  10.3  Hz 
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Figure  B-6.  Pressure  Field  Versus  Distance  at  15.1  Hz 
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Figure  B-8.  Pressure  Field  Versus  Distance  at  20.0  Hz 
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Figure  B-9.  Pressure  Field  Versus  Distance  at  22.4  Hz 
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Figure  B- 12.  Pressure  Field  Versus  Distance  at  29.7  Hz 
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Figure  B-13.  Pressure  Field  Versus  Distance  at  32.1  Hz 
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Figure  B-18.  Pressure  Field  Versus  Distance  at  44.2  Hz 
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Figure  B-19.  Pressure  Field  Versus  Distance  at  46.7  Hz 
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Figure  B-26.  Pressure  Field  Versus  Distance  at  63.6  Hz 
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Figure  B-27.  Pressure  Field  Versus  Distance  at  66.1  Hz 
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Figure  B-32.  Pressure  Field  Versus  Distance  at  78.2  Hz 
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